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ABSTRACT 

We analytically compute the effects that a pulsar's mass variation, whatever 
its physical origin may be, has on the standard Keplerian changes Arxep in the 
times of arrival of its pulses due to potential test particle companions, and on 
their orbital dynamics over long time scales. We apply our results to the plane- 
tary system of the PSR B1257-I-12 pulsar, located in the Galaxy at ~ 600 pc from 
us, to phenomenologically constrain a putative accretion of non-annihilating dark 
matter on the hosting neutron star. By comparing our prediction for Ar^^ ,j,j to 
the root-mean-square accuracy of the timing residuals 6{At) = 3.0 /is we find 
for the mass variation rate M/M < 1.3 x 10~^ yr~^. Actually, considerations 
related to the pulsar's lifetime, of the order of At ~ 0.8 Gyr, and to the currently 
accepted picture of the formation of its planets point toward a tighter constrain 
on the mass accretion rate, i.e. M/M < 10~^ yr~^. Otherwise, the planets would 
have formed at about 300 — 700 au from PSR B1257-I-12, i.e. too far with respect 
to the expected extension of 1 — 2 au of the part of the protoplanetary disk con- 
taining the solid constituents from which they likely originated. In fact, an even 
smaller upper limit, M/M < 10^^^ yr^^, would likely be more realistic to avoid 
certain technical inconsistencies with the quality of the fit of the timing data, 
performed by keeping the standard value M = IAMq fixed for the neutron star's 
mass. Anyway, the entire pulsar data set should be re-processed by explicitly 
modeling the mass variation rate and solving for it. Mo del- dependent theoretical 
predictions for the pulsar's mass accretion, in the framework of the mirror matter 
scenario, yield a mass increment rate of about 10~^^ yr~^ for a value of the den- 
sity of mirror matter pd^ as large as 10~^^ g cm~^ = 5.6 x 10^ GeV cm~^. Such 
a rate corresponds to a fractional mass variation of AM/M ~ 10~^ over the pul- 
sar's lifetime. It would imply a formation of a black hole from the accreted dark 
matter inner core for values of the dark matter particle's mass r/idm larger than 
3 X 10^ Gev, which are, thus, excluded since PSR B1257-I-12 is actually not such 
a kind of compact object. Instead, by assuming pdm ~ 10"^'' g cm~^ = 0.56 GeV 
cm~^, the mass accretion rate would be M/M ~ 10"^^ yr^^, with a fractional 
mass variation of the order of AM/M ~ 10~^^. It rules out mdm > 8 x 10^ Gev. 
Extreme values pdm = 1-8 x 10^^^ g cm~^ = 10^^ GeV cm~^ for non- annihilating 
dark matter in central spike may yield the constraint M/M < 10^^^ yr^^; over 
At = 0.8 Gyr, it rules out mdm > 12 Gev. 

Subject headings: gravitation— dark matter— planetary systems— pulsars: 
general— pulsars: individual, (PSR B1257+12)— extrasolar planets 
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Introduction 



1.1. Need for Dark Matter 



An increasing number of observations at galactic, extragalactic and cosmological scales, 
if interpreted in the framework of the presently accepted Newtonian/Einsteinian laws of 
gravitation, requires the existence of huge amounts of a peculiar kind of matter which 
does not emit electromagnetic radiation: the so-called Dark Matter (DM). It cannot be of 
baryonic nature. Indeed, measurements of the baryon density in the Universe using the 
Cosmic Microwave Background (CMB) spectrum and primordial nucleosynthesis constrain 
the baryon density to a value less than 5% of the critical density pcrit- Instead, the total 
density of clustered matter, obtained from Supernovas-based measurements of the recent 
expansion history of the Universe, CMB measurements of the degree of spatial flatness, and 
measurements of the amount of matter in galaxy structures obtained through big galaxy 
redshift surveys, is about 27% of the critical density. Thus, about 22% of it must exist in 
an exotic, unknown form. For reviews of b oth theor e tical a nd observation a l aspe cts of the 



DM paradigm, see, e.g., iBergstroml ( 120001 ): ICondold (J2004f ): iBertone et al.l ( 120051 ) 



1.2. Dark Matter accretion onto astrophysical objects 



If DM exists, it should be present in all astrophysical objects; it may both be there 
since their formation process and it may subsequently be accreted from the surrounding 
environment. In recent years much efforts ha ye been devoted to investigate the phenomenon 



Bertone Sz Fairbairn 2008: Kouvaris 2008: Sandin fc Ciarcelluti 20091: Ciarcelluti fc Sandin 



of possible capture of DM by neutron stars (IGoldman fc Nussinovlll989l: iGould et al.lll990 



2010 : Ide Lavallaz fc Fairbairnll2010l : iKouvaris fc Tinyakovll2010t IConzalez fc Reiseneggei 
20101) • Indeed, such compact objects should efficiently capture DM because of their high 
matter density. The content of DM depend on the nature of its particles, the type of 
hosting celestial bodies and their history. Moreover, new precise results from observations 
of neutron stars are becoming more frequently available. Thus, at least in principle, they 
are considered as potentially useful tools to independently constraining various aspects of 
DM models like density, cross section and mass of their particles. Such parameters are also 
crucial i n determining the c apabilities of several Earth-bas ed laboratory experiments like 
CDMSI (lAkerib et al.ll2003h . CDMSII (lAhmed et aDl2010h. D AMA/Nal (iBernabei et al. 



2003) and its successor DAMA /LIBR A ( IBernabei et al.l 120081 ). XENONIO (lAnde et al. 
2009[ i and ZEPLIN III (ISummerlliooJh aimed to directly detect DM. 
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1.2.1. Self- annihilating Dark Matter and some consequences of its accretion on 

astrophysical objects 

According to the widely popular Weakly Interacting Massive Particle (WIMP) scenario, 
DM annihilates with itself and interacts with the rest of the Standard Model (SM) via 
the weak interaction. The WIMP is typically defined as a stable, electri cally neutral, 
mass ive particle which arises naturally in supersymmetric SM extensions fJHaber fc Kane 
19851 ). A pair of WIMPs can annihilate, producing ordinary particles and gamma rays. 
Self-annihilating particles captured by neutron stars would contribute to alter their outward 
appearances because t he energy released in their a r inihilation wou l d affect their internal and 
surfa c e temperatures (IGoldrnan fc Nussinovlll989l: lKouvarisll2008l : Ide Lavallaz &: Fairbairn 
2010l : iKouvaris fc Tinyakovl I2OIOI : iGonzalez fc Reiseneggerl 120101 ). On the other hand, 
WIMPs do not steadily accrete onto neutron stars and do not substantially modify 
their inner structure by, e.g., inducing a massive DM core which may notably alter the 
gravitational collapse (see, instead, Section [1. 3. ip .. 



1.3. Non-annihilating Dark Matter. Mirror matter 

On the other hand, models of DM exist in w hich it does not undergo self-annihilation 
(JNussinovlll985l : lKaplarull992t iHooper et al.ll2005l ): it may happen, for example, if DM is 
made of fermions, without the corresponding antifermions, or if DM consists of bosons 
and carries one sign of an additive conserved quantum num ber, but not the opposite 
sign. Among such scenarios there is the mirror matter one flBlinnikov fc Khlopov 1982 : 



FootJl2008l : lBlinnikovll2010h . 



Khlopov et all 119911 : lKhlopovlll999l : lFootll2004al : IOkun1l2007l : 
The possible existence of su ch an exotic form of matter was e nvisaged for the first t ime 
in the pionee ristic works by iLee fc YangI (119561 ) and, later, bv iKobzarev et al.l (119661) and 
Pavsid (119741 ): the modern form of such an idea was laid out by lFoot et al.l (119911 ). Mirror 
matter arises if instead of (or in addition to) assuming a symmetry between bosons and 
fermions, i.e. supersymmetry, one assumes that nature is parity symmetric. In such a 
framework, in order to restore the parity symmetry violated by the weak interactions, the 
number of particles in the Standard Model is doubled in such a way that the Universe 
is divided into two sectors with opposite handedness that interact mainly by gravity. On 
the other hand, parity can also be spontaneously broken depending on the Higgs potential 
( Berezhiani fc Mohapatralll995l : iFoot et al.ll2000l ). While in the case of unbroken parity 



symmetry the masses of particles are the same as their mirror partners, in case of broken 
parity symmetry the mirror partners are lighter or heavier. In regard to the interaction 
mechanisms among such putative mirror particles, the forces among them are mediated 
by mirror bosons. Now, with the exception of the graviton, none of the known bosons 
can be identical to their mirror partners. Mirror matter can interact with ordinary matter 
in a non-gravitational way only through the so-called kinetic mixing of mirror bosons 
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with ordinary bosons or via the exchange of Holdom particles ( iHoldoml Il986a| ]bl) . These 
interactions can only be very weak. That is why r nirror particles have been suggested a s 



DM candidat es (JBlinnikov fc Khlopovlll982l . Il983l : iKolb et al.lll985l : iKhlopov et al.lll991 



Hodgeslll993[ ). Putting effectively constraints on the masses of such kind of stable particles 
in Earth-based accelerator experiments is mu ch more difficult than for self- annihilating 
candidates; in experiments like DAMA/Nal (JBernabei et al.l l2003l ) mirror DM would 
interact with ordinary matter via renormalizable photon-mirror photon kinetic mixing, 
leading to a recoil energy-dependent cross section. Actually, mi rror matter is one of the 
few DM candidates which can explain the positive DAMA/Nal (JBernabei et al.ll2003l ) dark 
matter signa l whilst still being consistent with the null results of other DM experiments 
( lFootll2004bl Jd). Another interesting feature of mirror matter is that it may overcome the 
difficulties that WIMP-based DM mo dels have in exp laining the opposite behaviors of DM 
in some colliding clusters of galaxies ( ISilagadzd l2009l ) . Indeed, while the behavior of Abell 



520 points towards a significant self-interaction cross-section for DM flMahdavi et al.l 120071) . 
the B ullet cluster (IE 0657-56), instead, behaves as a collisionless system ( IClowe et al. 
20061 ). As a result, the inferred estimate on the DM self-interaction cross section is well 



above the upper limit derived for the Bullet cluster and exceeds by many orders the cross 
section magnitude expected for WIMPs. On the contrary, mirror DM models exhibit a 
great er fiexibility; for ther n diverse behavio r of DM is a natural expectation ( ISilagadze 
20091 ). On the other hand, iBlinnikovl ( l2010l ) showed that the constraints on self-interaction 
cross-sections derived from observations of colliding clusters of galaxies are not real limits 
for individual particles if they form macroscopic bodies. 



1.3.1. Accretion of non- annihilating Dark Matter on astrophysical objects and some 

consequences of it 

Some consequences of accretion of non-annihilating DM onto r ieutron stars have 
more or less recently been inves t igated flBlinnikov fc Khlopovlll983t iBertone fc FairbairnI 



20081 : ISandin fc Ciarcellutill2009l : ICiarcelluti fc Sandinll2010l : Ide Lavallaz fc FairbairnI 120101 ). 



Basically, in this case a bulky mass of DM would accumulate inside a neutron star without 
self-annihilating; such kind of DM particles do not heat the star (cfr. with Section fl. 2. ip . 
The presence of such a DM core inside a neutron star may alter the usual mass-radius 
relation for su ch kind of astrophys i cal c ompact objects with potentially detectable 
consequences (ICiarcelluti &: Sandinll2010l ). Such an effect has been calculated in the 
framework of mirror DM, but it is qualitatively valid also for other kinds of non-annihilating 
DM that could form stable cores inside neutron stars. Moreover, if the mass accretion 
continues steadily, the growing inner DM bulk may reach its own Chandrasekhar mass 
l imit and collapse into a black hole, thus swallowing the hosting neutron star as well 
(Ide Lavallaz fc FairbairnI l2010l ). It must be noted that, in general, the Chandrasekhar 



"dark" mass M^^ would be smaller than in the usual case. Indeed, in terms of the Planck 
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mass Mpi (jde Lavallaz fc Fairbairnll201Cll ). 



M, 



(dm) 



Mh 



Ch 



m 



;i) 



dm 



where rudm is the mass of a DM particle which is larger than the ordinary nucleon mass 
"^nuc ~ 0.9 Gev in most models. Such a dramatic outcome, per se interesting, could, in 
principle, account for the unexplained gamma ray bursts observed in the Universe instead 
of resorting to the usual picture involving the coalescence of a neutron star with another 
compact object. 



1.4. Overview of the paper 

In this paper we will consider other consequences of the accr etion of non-annihilating 



DM by resorting to a specific scenario, i.e. the planetary system (IWolszczan fc Frail! Il992 



Rasio et al.lll992h hosted by the PSR B1257+12 pulsar flWolszczanlll990l ) situated in the 
Galaxy at high Galactic latitude at less than 1 kpc from us. In particular, we will look at 
some effects connected with the planets' orbital dynamics induced by a putative increment 
of non- annihilating DM experienced by the pulsar to derive bounds on such a phenomenon 
which are, then, used to constrain the resulting dark core. A somewhat analogous study was 
performed with stellar motions around the Galactic Center to constrain DM annihilation 
propo sed to explain the Te V gamma rays emanating from such a regio n of the Mi l ky 



Way (JHall fc Gondoldl2006l ). Concerning mirror matter and exoplanets, iFootI ( 1l999l . l200ll ) 
suggested that the several close-in planetary companion s of main-sequence stars discovered 
so far may be made up primarily of such a form of DM. iFoot fc Silagadzd (J200l[ ) looked for 
mirror planets in our solar system itself. Our approach can, in principle, be extended also 
to other similar scenarios involving a pulsar harboring compact or standard companion(s). 
In Section |2] we will put phenomenological, model-independent constraints on M/M from 
the analysis of the impact that it may have on the standard Keplerian variation of the times 
of arrival (TOAs) of pulsar's pulses (Section 12. ip by comparing the present-day accuracy in 
their root-mean-square residuals with the analytically computed variation of TOAs due to 
M/M (Section 12. 2p . We will also use some considerations on the age of the pulsar's system 
(Section 12. 3p . Then, we will consider the global reduction of the spatial extension of the 
orbits of the planets of PSR B1257+12 due to its supposed mass accretion by contrasting 
our predictions of such an orbital shrinkage with the currently accepted picture of the 
birth of the system (Section 12. 4p . In Section |3] we will, first, compare our constraints to 
some recent predictions for the pulsar's mass accretion in terms of the mirror DM scenario 
(Section [3. ip . Then, we will trace some consequences in terms of dark core collapse (Section 
13. 2p . Section H] is devoted to the conclusions. 
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2. Planetary orbital effects of the mass accretion onto PSR B1257+12 



An ideal laboratory to study certain consequences of a putative accretion of non- 
annihilating DM on a neutron star in a purely phenomenological and model-independent 
way is represented, in principle, by a system hosting a pulsar orbited by one or more 
companions of planetary size; PSR B1257+12 a nd its three small Earth-siz ed planets 
represent one of such scenarios alo 600 ± 100 pc ( Konacki fc Wolszczanll2003l ) from us, with 
Galactic longitude and latitude / = 311 deg, b = 75 deg, respectively. Indeed, in such cases 
one can look at various effects induced by the pulsar's mass variation on the orbital motions 
of its planets. Actually, the strategy devised below applies to any kind of putative mass 
variations of the primary. 



2.1. The Keplerian change in the times of arrival of the pulsar 

The direct observable is the change At in the pulsar's TOAs due to its orbiting 
partners. The usual periodic variation Arxep of TOAs resulting from the p urely Keplerian 



motion of a pulsar around the center of mass of the star-planet system is (JKonacki et al. 
2000h 



A 



TKep 



X 



(cos E — e) sin uj + a/1 — e^ sin E cos uj 



(2) 



where x = a sin I /c, [x\ = T, is the projection of the pulsar's semimajor axis a with respect 
to the center of mass, e is the eccentricity of the orbit, E is the eccentric anomaly and cu 
is the argument of periastron. Here c denotes the speed of light in vacuum, [c] = L T~^. 
The inclination / is the angle between the orbital angular momentum and the line-of-sight 
from the pulsar to us. The semi-major axis a characterizes the size of a Keplerian ellipse: 
dimensionally, [a] = L. The eccentricity e is an adimensional parameter which fixes the 
shape of a Keplerian ellipse. It is a non-negative real number which can assume all values 
within < e < I, where e = corresponds to a circle. The eccentric anomaly E can be 
regarded as a parametrization of the polar angle in the orbital plane. The longitude of 
periastron uj is an angle in the orbital plane which determines the position of the point of 
closest approach, generally dubbed periapsis, with respect to a reference direction which 
is customarily assumed coincident with the line of the nodes. The line of the nodes is 
the intersection between t he orbital plane and the plane of the sky, chosen in this case as 
reference plane. See, e.g., ( lRoyll2005l ) for basic concepts on orbital mechanics. 



For a better comprehension, let us point out that eq. ([2]) yields the instantaneous value 
of ArKep, corresponding to a given value of E: indeed, there is a mapping between the 



^Such 



an 



estimate for the distance is based 



on 



tic electron distribution model by iTaylor fc Corded ( 119931 ) 



the Galac- 
See also 



http://vizier.u-strasbg.fr/viz-bin/VizieR-S7PSR%20B1257%2bl2 



eccentric anomaly and time given by 



E = M + J2(lj M^e) sm{iM), M = n{t- tp), 



(3) 



where Ji{ie),i = 1, 2, . . . are the Bessel functions of first kind (1Watsonlll966l ). Ai is the 
mean anomaly ( JRoyll2005l ). which is a parametrization of time, and tp is the time of passage 
at periastron. The series of eq. ([3]) co nverges for all values of e < 1 like a georn etric series of 
ratio (eexp \/l — e^) / (l + Vl — e^) ( IWintnerlll94ll : Ida Silva Fernandeslll994l ). During an 



orbital revolution E spans an angular interval of 27r, in such a way that the timing r does 
not remain constant, as it happened if the pulsar was not perturbed by its companions, 
but exhibits a time-dependent, harmonic variation At{E) which reveals the existence of 
other bodies in the s ystem. When a pulsar has A^ companions, the TOA variations become 
flKonackiet alJl2000[ ) 



N 



ArKep = J2 






l^cos Ej 



sm Uj + 



1 — e^ sin Ej cos Uj 



(4) 



i=i 



2.2. The non-Keplerian perturbations of the change in the times of arrival of 

the pulsar due to its mass variation 

The further change in TOAs with respect to the purely Keplerian one of eq. ([2]) due 
to a generic orbital perturbation caused by a small force which deviates from the largest 
Newtonian two-body, pointlike monopole acceleration GM/r"^ causing the well-known 
Keplerian motion can straightforwardly be worked out from eq. ([2]) itself by differentiating 
it with respect to the orbital parameters which undergo slow time-variations induced by 
the perturba tion conside red. In the case of the PSR B1257+12 system it has recently 
been shown (jloridl2010bl ) that standard non-Keplerian dynamical effects (departures from 
sphericity of the pulsar, IPN, Schwarzschild-like corrections of order 0(c~^)) which may 
cause departures from the main Keplerian picture are negligible, given the present-day 
accuracy in measuring some orbital characteristics like the orbital periods Ph of the planets. 
To better understand such points it is useful to note that the ratios of the pulsar's radius 
{R ~ 10 km) to the planetary orbital separations (about 0.2 — 0.5 au, see Tabled]) are of 
the order of just 3.5 — 1.5 x 10~^. Thus, the point mass approximation is well justified, 
and the corrections induced by general relativitMj to th e values of both the mass and 
current multipole moments (JLaarakkers fc PoissonI Il999l ) of the external gravitational 



^Since neutron stars are compact objects with strong internal gravity, their gravitational 
fields must be fully described within the framework of general relativity by using appropriate 
equations of state for their ultradense matter. 
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field of the neutron star are completely neKligib le as far as so distant orbital motions are 
concerned (jSibgatullinI l2002l : iBejger et al.ll2010l ). The TOAs are, in principle, also affected 
by effects concerning the propagation of the pulsar' s electromag netic waves in the distorted 
space-time (gravitational red-shift, S hapiro delav) fIStairsI 120031) usually accounted for by 
the post-Keplerian parameterqj 7, r ( Damour fc Deruellelll986l ). However, inserting the 
figures of Table [1] into the analytical expressions for 7 and r shows that they are completely 
negligible in the present case with respect to the accuracy in measuring the TOAs (see eq. 
( ITOll later). Thus, we will not further consider them in the following. 

In the specific case of a putative pulsar's mass variatioiiu with percent rate M/M the 
perturbation to Arxep is 



A^ 



M/M 



\ dx 



^■^ ^ Ax + r^^Ue ■ ^^^^^'^P 



V de 



V dE 



--(^l-' (^) 



^The parameter 7 is the amplitude of the combined effect of the special relativistic time 
dilation and the gravitational red-shift, while r is the amplitude of the Shapiro propagation 
delay caused by the gravitational field of the pulsar's companion. For an insightful qualitative 
description of such phenomei ia, which helps in understanding why they are negligible in the 
PSR B1257 + 12 system, see iKrameJ fl2010h . 



^As recently shown both analytically and numerically by llorid (J2010al ). such an effect 
is non-Keplerian in the sense that the resulting orbital motion of a test particle about the 
mass-varying primary is not a closed Keplerian ellipse. Note tha t this is a cl assical orbital 
effect, being the general relativistic orbital one totally negligible ( llorid l2010a[ l. 



Table 1: Relevant orbital param e ters o f the three planets A,B,C, of PSR B1257+12, from 
Table 2 of iKonacki fc Wolszczanl (J2003[ ). Here ap are the planet semimajor axes. Figures in 
parentheses are the formal 1 — a uncertainties in the last digits quoted. The standard value 
M = IAMq for the pulsar's mass has been kept fixed in deriving them, i.e. it has not been 
included in the list of the parameters to be solved-for in the fit of the pulsar's timing data. 





A 


B 


C 


X (ms) 


0.0030(1) 


1.3106(1) 


1.4134(2) 


ap (au) 


0.19 


0.36 


0.46 


e 


0.0 


0.0186(2) 


0.0252(2) 


n(d) 


25.262(3) 


66.5419(1) 


98.2114(2) 


uj (deg) 


0.0 


250.4(6) 


108.3(5) 
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with 



( dATKe 



dx 
de 
dE 

9ArKep 

duj 



yjl — e^ sin E cos u + (cos E — e) sin u, 



-x|^^i£££^ + sina; 



X 



X 



(Vl — e^ cos E cos w — sinii^sino;) , 



(6) 



[(cos E — e) cos a; — Vl — e^ sin i? sin w] . 



In eq. (|5]) the variations of a, e, -E, u d ue to M/M appear; J is left unaffected by the mass 
variation of the primary ( lloriol 2010al): thus . Ax = Aa sin J/c. They have been exphcitly 
computed as functions of E by lloriol ( l2010af l: 



Aa 



m\ 2ae / sin, 
M I n \ 1 



E-E cos E 
-ecos E 



) 



A^ = -(S)^(^ 



E-E COS E 
- e cos B 



)^ 



(7) 



AE 



Aw 



f ) i [^(i?) + S(E) + C(E)] 



M 



{l+e){l-cosE)-EsmE 
1— ecos E 



where the coefficients A, B, C of the variation of the eccentric anomaly are 

A(E) 



E'^+2e{cosE-l) 
1— ecosB ' 



I B{E) 



C{E) 



l-e^ 



[(l+e)(l-cos-E)-Esin£;] 
(l-ecosE)2 ■ 



(1— e^) sin E {sin E— ecos E) 
[l—ecosE)'^ 



In eq. ^ n = 27r/Pb is the Keplerian mean motion and Pb is the Keplerian orbital period. 
For those readers not specifically acquainted with the methods of celestial mechanics it may 
be useful to point out some features of eq. ([7]). They represent the changes in a,e,E,u 
induced by the considered perturbation at a given instant, to which correspond a given 
value of the eccentric anomaly E. Thus, they are not to be intended as the variations per 
orbit. Thus, eq. ([5]) represents the instantaneous change in TOA, so that its variation per 
orbit must be evaluated by taking the difference between Atj^j^j computed at -E = 27r 
and ATjyjn.]: computed at E = 0. Incidentally, note also that the right-hand-sides of eq. 
([7]) have correctly the same dimensions of the left-hand-sides. Indeed, [n~^] = T, while 
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[M/M] = T\ After k = l,2,3,... revolutions, the TOA variation A^^ due to M/M is 



^£/M = ^^a//m(2^^) - A^A//A/(0) = M ^ 1 ^bX 



M 
M 



1 + e 

k2TT\l cos a; — (1 — e) sinw 

1 — e 



■ (9) 



As stated before, a useful application of the previously obtained results is represented 
by the planetary system of the 6.2 ms PSR B1257+12 pulsar. Such a neutron star was 
discovered in 1990 during a high Galactic latitude s earch for millisecond pulsars with the 
Arecibo radiotelescope at 430 Hz (jWolszczanI Il990l ) . Two years later, PSR B1257+12 
turned out to be orbited by at least tw o Earth-sized planets-dubbed B and C-along almost 
circular paths ( IWolszczan &: Fraillll992l ). In 1994 it was announced the di scovery of a third. 
Moon-sized pl anet-named A - in an inner, circular orbit (JWolszczanI 19941). Its pre s ence, 
questioned bv IScherer et al.l (Il997l ). was subsequently confirmed in iKonacki et al.l (Il999[ ): 
Wolszczan et all fl2000bf )" 



Concerning the accuracy in timing PSR B1257+12, covering 12 yr, a detailed 
de scription of the data acq uisition and the TOA measurement process can be found 



Wolszczan et al.l ( l2000al ). The final post-fit resid uals for daily-averaged TOA s are 



m IWolszczan et al.l l izuuual ). ine nnai post-nt resia uais tor dany-averaged lUA 
characterized by a root-mean-square (rms) noise ojj (IKonacki fc Wolszczanll2003l ) 



(5(Ar) = 3.0 /iS. 



(10) 



By equating eq. ([9]) for the three planets to eq. (fTOj) it is possible to obtain upper 
bounds for a putative mass variation experienced by PSR B1257+12; the results are in Table 
El and do not depend on any specific model of non-annihilating DM. Anyway, it should 
be pointed out that they should be considered just as preliminary, order-of- magnitude 
evaluations; actually, the entire pulsar timing data set should be re-processed by explicitly 
modeling the pulsar's mass variation as well, and a dedicated solve-for parameter should be 
estimated in a least-square sense. From Table |2] it can be noted that the tightest bound 
amounts to 

< 1.3 X 10"" yr-\ (11^ 



M 



In Section 12.31 and Section 12.41 we will independently check eq. ( ITTj) by inspecting its 
compatibility with other effects connected with the dynamical history of the PSR B1257+12 
system. 



^ The relativistic contributions of order 0{c ^) to At due to gravitational time delay and 
the Shapiro delay are far smaller than eq. (TTOj) . 
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2.3. Considerations from the age of the pulsar's planetary system 



Concerning the planets of PSR B1257+12, it is just the case to review here that 

such a peculiar planetary system got formed after the death of the star progenitor of PSR 

B1257+12, not at it s birth as a main sequence star, as in the usual planetary formation 

(jPodsiadlowskil Il993l ) . Thus, the age of PSR B1257+12 as a neutron star cannot be smaller 

than ~ 10^ yr, which is roughly the time required to form planets around millisecond 

pulsars from pro toplanetary disks starti ng from the epoch of the birth of the pulsar itself 

as stellar corpse (ICurrie fc Hansenll2007l ): this implies that, by setting At = 10^ yr, eq. (fTTj) 

yields 

AM ^ , , 

Actually, such a constrain seems to be too large. Indeed, apart from the fact that it 
would imply that, in the framework of the non-annihilating DM accretion scenario, PSR 
B1257+12 should have already become a black hole, at least for certain values of DM 
particle's mass (see Figure 8 of Ide Lavallaz fc Fairbairru (J2010[ ) and eq. (1271) below), on the 
other hand its mass should now be as large as l.SSM© with respect to an assumed standard 
value of M = \AMq at its birth. But the post-fit residuals in the TOAs of PSR B1257+12 
are statistically compatible with zero, i.e. they do not exhibit anything anomalous at a 
statistically significant level. To this aim, it should be noted that such residuals have 
been obtained just for M = IAMq. Stated differently, if the mass of PSR B1257+12 was 
really larger than the standard value by roughly 13%, the timing residuals, constructed 
just by keeping fixed the pulsar's mass to the standard value, should have retained some 
statistically significant non-zero features, which is not the case. Thus, we conclude that 
\M/M\ should be actually smaller than 10~^ yr~^ by likely one-two orders of magnitude. 
In fact, a more quantitatively precise evaluation of it would require a re-processing of the 
pulsar timing data with different values of its mass to check the level at which departures 
from the standard value M = IAMq cease to affect the TOAs post-fit residuals. It must 
be noted that such considera tions may be quite optimistic since the age of PSR B1257+12 
may be as large as ~ 1 Gyr ( Brvden et al.l 120061). The age of the pulsar is estimated from 
its spin-down timescale V/2V (JBryden et al.ll2006l ) to be < 1 Gyr since the pulsar's spin 
period is P = 6.2 ms and its variation is "P = 1.1 x 10^^^ ms s~^ ( jKonacki &: Wolszczan 



Table 2: Order-of-magnitude upper bounds for the mass accretion rate of PSR B1257+12 
from eq. ([9]) applied to the three p lanets of the system. An ac curacy of 6{At) = 3.0 /^s in 
the TOAs residuals has been used (JKonacki fc Wolszczan! |2003| ). A timing span of At = 12 
yr has been assumed corresponding to approximately /ca ~ 170, fee ~ 66, /cc ~ 45 orbital 
revolutions for A,B,C. They have been obtained in a phenomeno logical, model-independent 
way. 



A 



B 



C 



\M/M\ (yr-i) 7.6 x 10"^ 1.3 x 10"^ 1.9 x 10" 
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20031 ). In this case, by repeating the previous calculations, M jM should be, perhaps, of 



the order of 10~ yr~ or less to avoid striking contradictions with the lacking of likely 
statistically significant patterns in the TOAs residuals obtained for M = 1.4Mq. All in all, 
it should be recalled that to be more quantitative the phenomenon of mass variation of 
the pulsar should be explicitly modeled in the dynamical force models used in the timing 
processor systems and explicitly solved-for in the consequent re-analysis of the timing data. 



2.4. The effect of the pulsar's mass variation on the orbits of its planets 

throughout their history 



An independent check of the previous considerations can be obtained by inspecting 
another consequence of the mass variation of the primary on the orbital motions of 
its test particle companions: the change in the size of their orbits. Indeed, if PSR 
B1257-I-12 was really accreting its mass since its birth, the spatial extension of its planetary 
system should have been larger than now when it formed. Although no observationally 
determined quantities related to such a putative variation of its size are available for the 
PSR B1257+12 system, it is interesting to inspect the consequences that such an orbital 
shrinking may have on timescales comparable to the pulsar's lifetime a nd compare them to 



the currently accepted picture of the formation of its planetary system fiBryden et al.l 12006 



Currie fc HansenI 120071 ). From the expression of the Keplerian planet's astrocentric distance 
(JR,ovll2005h 



it follows 



a(l — ecosE"), 



Ar{E) = (1 - e cos E) Aa - a cos E Ae + ae sin E AE. 



(13) 



(14) 



It represents the instantaneous departure of the astrocentric distance with respect to the 
unpertu r bed, K eplerian one. Note that eq. (JHj) agrees with the results obtai ned in, e.g., 
Casottd (119931 ). In the case of a mass variation, from eq. (I71)-eq. (jH]) we have (jloridl2010al l 



Ar,,/,,iE)={^]^[ViE)+TiE)], 



(15) 



with 



V{E) 
HE) 



-2(sin E-EC0SE) + -'"^[^;+^-(-^"l)] _ (l-e^)sin^i.(sinZ; ecosE) 
V / 1— ecosB (1— ecosB)-^ 



( i-ecosi? ) {cosE{smE- E COS E) + smE 



(l+e) (l-cos E)-E sin E 
1— ecosB 



(16) 
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fc orbital revolutions is 



It turns out from eq. ( IT5l) and eq. ( IT6l) that the shift in the star-planet distance (JJ^v ,„ after 



4?/M = ^^M/M(2A:7r) - Ar,,/,,(0) = "M ]^ 1 ^ba(l - e). (17) 

As expected, eq. flT7|) tells us that the orbit gets smaller for a mass increase of the primary, 
i.e. 5[- , < for M /M > 0; note that here a denotes the relative star-planet semimajor 

axis whose values are listed in Table [H The application of the upper bound of eq. (fTTl) to 
eq. (fT7|) over a past time span of the order of the system's lifetime, i.e At ~ 1 Gyr, yields 
the implausible results listed in Table O A,B,C, should have formed at hundreds au far 
from the pulsar to have reached nowadays their present astrocentric distances under the 
action of a mass accretion experienced by their primary as large as that of eq. fllip . A 
mass accretion rate three orders of magnitude smaller than eq. fITT]) for PSR B1257-I-12, i.e. 
of the order of 10~^ yr~^, would, instead, yield primeval astrocentric distances as large as 
about 0.7 au, which is substantially in agreement with the maximum extension of ~ 1 — 2 
au of the p art of the protoplanetar y disk containing solid materials from which A,B,C likely 
originated fICurrie fc HansenI 120071 ). On the other hand, an age of the pulsar as large as 



about 1 Gyr would yield AM/M ~ 0.8 — 1, so that the same problems encountered in 
Section 12.31 would occur. Also in this case, a mass accretion rate two orders of magnitude 
smaller, i.e. M/M ~ 10^^^ yr^^ which would be well compatible with the distances at 
which the planets should have formed, would cure them. 



3. Confrontation with some theoretical scenarios 

3.1. Predictions of some scenarios for mirror matter accretion 

Let us, now, consider the predictions for the accretion of non-annihilating DM. It is 
just the case to briefly recall that while WIMPs may heat a neutron star without creating 
a particularly massive dark core inside, on the contrary, mirror particles do not heat a 
neutron star and may form a dark core with substantial mass that modifies the structure 
of the star. Reasoning in terms of mirror matter, its distribution in galaxies is expected 



(k) 

Table 3: Variations S)-J,,r of the astrocentric distances, in au, of the three planets of PSR 

B1257-M2 due to M/M = 1.3 x 10"^ yr"^ A time span of At = -1 Gyr has been assumed 
corresponding to approximately /ca ~ 1 x lO^'', /cb ~ 5 x 10^, /cc ~ 3 x 10^ orbital revolutions 

for A,B,C. 

ABC 

Ig^^(au) -304 -565 -717 
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to be non-homogeneous since it should form complex structures similarly to as ordinary 
baryons do. Thus, the accretion rate of mirror matter by neutron stars should depend on 
the location and history of each star. Of course, also the s t ructu re of the hidden mirror 
sector, which is unknown, is relevant. ISandin &: Ciarcellutil ( l2009l ) yield 

M = 10^ kg s"^ (18) 

as possible u pper bound of accretion of mirror matter onto a neutron star. Such an estimate 
is based on (IShapiro &: Teukolskylll983l ) 



M 



Pdn 



10 km s 



M 



kg s' 



(19) 



10-24 g cm- VV ^dm J \Mq^ 

which has been derived in the hypothesis that the distribution of mirror particles is isotropic 
and monoenergetic at large distance from the neutron star. In obtaining eq. f lTS]) from eq. 
f in?]) it has been assumed that the density of the interstellar medium in the mirror sector 
may be up to one order of magnitude larger than the value / 3Hm = 10-^^ g cm-^ ~.^-6 x 10^ 
Gev cm-"^ of the ordinary baryonic giant molecular clouds ( ISandin fc Ciarcellutil l2009l ). 
who tacitly use fdm = 10 km s~^ and M = Mq as well. Note that such a figure is orders 
of magnitude larger than the currently accepted value of the dark halo density at about 
10 kp c from the Galac t ic Center which is of the order op ^ 10-^*^ g cm" ^ = 0.56 Gev 
cm"^ (JBergstroml I2OOOI : iBertone et al.ll2005l: Ide Lavallaz fc Fairbairnll2010l ). Incidentally, 
let us note from Figure 1 of Ide Lavallaz &: FairbairnI ( 20101 ) that all the Einasto DM 
density profiles substantially converge to such a figure for PHm at Galactocentric d istances 



larger than 10 kpc, so that the uncertainty of 100 pc (JKonacki fc Wolszczanll2003l ) in the 
heliocentric distance of PSR B 1257+ 12 is of no concern for us. It is just the case of 
stressing that the estimate from eq. ( IT9l) has to be intended as referred to a local quantity. 
Indeed, the average DM density in the mirror matter scenario is comparable to that of other 
DM models. Contrary to them, mirror matter is believed to have the possibility of forming 
structures having higher density than the average DM one. In the case of PSR B1257+12, 
by assuming that it steadily accreted DM throughout its lifetime 

At = V/2V = 0.893 Gyr, (20) 

eq. (ITSi) would translate into a mass accretion rate of 



M 

M 



1.1 X 10"^^ yr"\ 



(21) 



It is compatible with the phenomenological constraints previously obtained for the PSR 
B1257+12 system. According to eq. (EI]) and eq. (EOj), PSR B1257+12 would have accreted 
a mass fraction 



AM 



1.0 X 10" 



(22) 



g cm ^ corresponds to 5.6 x 10^^ Gev cm ^. 
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during its lifetime. If, instead, we use eq. (fT9|) with pdn 
it yields an accretion rate of approximately 



10 ^^ g cm ^ 



0.56 Gev cm ^, 



M 
M 



2 X 10"^^ yr"^ 



corresponding to a fractional mass variation of PSR B1257 +12 over its lifetime of 

AM 



M 



1.9 X 10 



-14 



(23) 



(24) 



Also in this case, such figures are comp atible with the constraint s phen omenologically 
obtained. Let us mention that, recently, iGonzalez fc Reiseneggerl ( l2010l ) yielded a mass 
accretion rate for pulsars in the solar neighborhood of the order of 



M < 9 X 10"^^ Mq yt-\ 



(25) 



although it seems that they considered self-annihilating DM. Let us note that that a mass 
accretion rate of about 10~^^ yr~^, which would be in agreement with both the planetary 
formation scenario and the fit of the timing data for PSR B1257+12 (Section l2.3J -Section 
12. 4p . could be obtained from eq. ( TT9l) by usi ng the extreme limit Pd m = 1-8 x 10"^'^ g 
cm~^ = 10^^ GeV cm^^ of the predictions by lBertone fc MerritI ( 120051 ) for non-annihilating 
DM in a central spike. 

Another approach that, in principle, one could follow consists of performing specific 
numerical simulations coming up with a quantitative estimate of the probability that 
PSR B1257+12 eventually emerged where it is now located starting from some detailed 
evolutionary models of the distribution of the mirror matter in the Galaxy taking into 
account issues like its form, dumpiness, density of dumpiness, etc, and by using different 
values of its density. This is beyond the scope of the present paper. Let us mention that 
such relatively "fine-graining" simulations of mirror matter distribution at the level of 
the Milky Way sub-structure have not (yet?) been impleme nted; some work has been 
performed at a more general level concerning gal actic haloes ([Mohapatra fc Teplitzll2000f). 
and cosmo logical large scale structure formation (jlgnatiev fc VolkasI 120031 ; iBerezhiani et al 



20051). 



3.2. Consequences in terms of dark core collapse for WIMPs only 

Let us interpret our results in terms of dark core collapse. In such a framework, the 
limiting case of the mass variation of a neutron star occurs when its putative dark core has 
reached its Chandrasekhar mass, i.e. 



AM^™^^) . [M^cr^ + M]-M M, 



(dm) 
Ch 



M 



M 



M 






1.1655 Gev^ 



m 



(26) 



dm 
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According to the values for m^^ reported in Table HI the maximum and minimum values 
for eq. ( l26l) are, by assuming M = IAMq, 



AM M') 
M 



1.1 X 10 



-16 



AM (coll) 
M 



(27) 



= 4.6 X 10"^. 
Since, evidently, PSR B1257+12 has not (yet?) become a black hole, it must be 




At < 






(28) 



non— ann 



According to our model-dependent estimates on AM/M of eq. fl2^ and eq. f l2^ . the 
condition of eq. fl25]) is clearly satisfied for the smallest admissible values of mdm connected 
to the maximum figure in eq. (!27|) . To this aim, let us incidentally note that in the mirror 
matter scenario the lowe r limit for ? 7 irim ca n actually be much smaller than that quoted 
here in Table HI see, e.g.. lFoot et al.l ( l2000l ). On the contrary, eq. ( 128!) does not hold for the 
largest supposed values of mdm yielding the smallest values in eq. fl^ . Thus, eq. ( 125]) is 
able to put upper bounds on mdm- To this aim, eq. ( ETj) allows to rule out DM candidates 
with masses larger than 



mdr, 



M|, 



\ M IM/M 



3 X 10^ Gev 



At 



(29) 



non— ann 



because they would yield limiting core-collapse fractional mass variations smaller than 
eq. (12^ . Instead, the much smaller figure of eq. (1231) yields an upper bound on the DM 
particle's mass 

mdm = 8 X 10^ Gev. (30) 

Moving to the phenomenological, model-independent constraints on M/M of Section 
12.21 Section [231 it is apparent that the y do not fulfil the conditio r i of eg . (125|) . at least for the 
values of mdm considered in Figure 8 of Ide Lavallaz fc FairbairnI (120101 ). However, as already 



Table 4: Values, in Gev, of some mass energies used in the text. Mq is the mass of the Sun, 



(min) (max) 



Mpi is the Planck mass, "n^dm ' '^dm ^^^ ^^^ r ninimum and maximum masse s of a non- 
annihilating DM particle according to Figure 8 of Ide Lavallaz fc FairbairnI ( 20101 ). Actually, 



in the mirror matter scenario rn^ can be much smaller than 50 Gev ( Foot et al.l 120001 ). 



Mq (Gev) Mpi (Gev) 



m 



dm 



(Gev) 



m 



(max) 
dm 



(Gev) 



1.1157 X 10^7 1.22105 X 10^^ 50 



10« 



-18- 



noted, i n the mirror matt er scenario the mass of the DM particles can be much smaller than 



50 Gev (JFoot et al.l 120001 ) . A mass accretion rate of 10^^ yr~^, which is compatible with the 
planetary formation history (Section I2.4p . would yield over At = 0.8 Gyr an upper bound 
on the DM particle's mass of 1.2 GeV. If, instead, we take M jM = 10~^^ yr~^, which would 
yield just a 0.8% departure of the pulsar's mass from its standard value over At = 0.8 Gyr 
(Section l2.3J -Section 12. 4p . the upper bound is rridm < 12 Gev. 



4. Summary and conclusions 

A neutron star in the Galaxy may accrete non-annihilating DM at a rate up to about 
M = 10^ kg s~^ for pdm ~ 10"^^ g cm"'^ ~ 5.6 x 10^ GeV cm~^ according to some authors; 
if, instead, we assume pdm ~ 10^^^ g cm"'^ ~ 0.56 GeV cm""^, the mass accretion rate 
would be about Af = 1 kg s"^; extreme values like pdm ~ 10"^'^ g cm~^ ~ 10^^ GeV cm~^ 
are also possible according to other researchers, yielding a mass accretion rate of the order 
of 10~^^ yr~^ . The steady accumulation of non- annihilating DM inside a neutron star 
may yield the formation of an inner dark core. If it reaches its own Chandrasekhar mass 
M^j^™ ~ Mpj/m^jjj, which should be different from that made of ordinary baryons in view of 
the expected different mass mdm of the the DM particles with respect to standard nucleons, 
such a core may collapse into a black hole, thus destroying the hosting neutron star. We 
used the PSR B 1257+ 12 millisecond pulsar, located at ~ 600 pc from us and hosting a 
planetary system of three Earth-sized companions, to put phenomenologically constraints 
on such a putative mass accretion rate M/M by looking at some dynamical orbital effects 
affected by such a phenomenon. We also exploited the expected lifetime of the neutron 
star. 

In the case of the PSR B1257+12, the aforementioned predicted mass variation 
rates would correspond, for M = IAMq, to M/M =~ 10"^^ yr'^ M/M ~ lO'^^ yr-\ 
and M/M ~ 10~^^ yr"^, respectively. They are all compatible with the upper bound 
M/M < 1.3 X 10~^ yr~^ phenomenologically obtained by comparing the analytically 
calculated perturbations by M/M on the standard Keplerian variations Arxep of the times 
of arrival of the pulsar's pulses due to the presence of its planets to the root-mean-square 
residuals 6{At) = 3.0 ps of the TOAs. Smaller phenomenological constraints by some 
orders of magnitude {M/M ~ 10~^ yr^^) come from the confrontation of the predicted 
global shrinking of the orbits of the planets during the pulsar's lifetime At ~ 0.8 Gyr to 
the currently accepted picture of their formation from a protoplanetary disk whose part 
containing solid particles extended just for 1—2 au; for M/M ~ 10~^ yr~^ the planets would 
have formed at approximately 300 — 700 au from the pulsar. Given the pulsar's lifetime, in 
order to avoid possible contradictions with the fact that the standard value of the pulsar's 
mass M = IAMq, kept fixed in its timing data fitting, did not destroy the goodness of 
the fit at a statistically significant level, a smaller mass variation rate M/M < 10^^^ 
yj,-i j^g^y i^g considered more plausible. However, to be more quantitative, the entire 
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pulsar timing data set should be re-processed by explicitly including the effect of M jM in 
the dynamical models used. According to some predictions of the non-annihilating DM 
scenario, during its lifetime PSR B1257-I-12 should have accreted a mass fraction as large 
as AM/M = 1.0 X 10"^ (pdm ~ 10"^^ g cm'^ ~ 10*^ GeV cm'^), or AM/M = 1.9 x 10"^^ 
(pdm ~ lO"^'' g cm~^ = 0.56 GeV cm~'^), respectively. Since it is not (yet?) a black hole, 
such figures exclude values for the DM particle's mass larger than 3 x 10^ Gev and 8 x 10^ 
Gev, respectively. Instead, M jM ~ lO"^'^ yr~^ and the pulsar's lifetime would imply a 
fractional mass accretion AM/M which rules out m^-c^ > 12 GeV. 

The approach followed here may be, in principle, extended to other scenarios involving 
one pulsar hosting compact or planetary companions. 
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